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PREFACE 


~  7/  '/  C  J 

This  technical  note  provides  a  current  description  of  U5AFETAC  Soil  Moisture  and 
Agronet  programs.  Changes  in  these  programs  over  the  past  two  years  have  left  the 
previous  report,  USAFETAC  TN  77-3,  outdated  in  many  of  its  particulars.  While  I  have 
relied  on  TN  77-3,  1  have  not  merely  patched  up  the  previous  report.  The  material 
has  been  updated,  expanded,  and  reorganized.  For  each  program,  user  products  are 
described  first  and  followed  by  sections  on  the  general  methods,  the  calculations, 
and  lastly,  the  production  programs .  This  way,  the  reader  does  not  have  to  progress 
through  the  details  of  how  the  programs  work  to  find  out  what  they  do.  But  I  hope 
someone  wi 1 1  want  to. 
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INTRODUCTION 


The  U5AF  Envi ronnental  Technical  Applications  Center  (USAFETAC)  prepares 
neteorol ogi ca  1  analyses  tailored  for  use  by  government  agencies  to  monitor  crop 
growing  conditions  over  much  of  eastern  Europe,  Asia,  and  the  United  States. 
USAFETAC's  support  to  these  customers  is  provided  under  two  projects,  known  as  "Soil 
Moisture"  and  "Agronet." 

The  Soil  Moisture  program  has  been  active  since  1958.  It  uses  C.  U. 
Thornthwaite 's  method  (Thornthwaite  and  Mather,  1957)  to  compute  water  balance. 
Thornthwaite 's  method  requires  only  precipitation  and  temperature  data  as  routine 
i nput --el enents  readily  available  from  surface  observations.  The  approach  is  simple, 
and  it  has  been  providing  useful  information  for  over  20  years. 

The  Agronet  program  became  operational  in  1974.  It  was  developed  as  interest  in 
worldwide  crop  production  grew  in  the  early  1970s  and  data  became  available  at  the 
Air  Force  Global  Weather  Central  (AFGWC)  to  meet  the  requirements  of  a  more  complex 
model .  Among  its  advances  over  the  Soil  Moisture  program,  the  Agromet  program  uses 
AFGWC 1  s  Automated  Cloud  Analysis  Model  (3DNEPH)  to  estimate  the  effect  of  clouds  on 
the  radiation  balance  and  uses  Penman's  (1956)  evaporation  formula. 

The  models,  procedures,  and  products  are  described  in  Chapters  1  (Soil  Moisture) 
and  2  (Agromet).  Some  basic  terms  are  defined  in  the  Glossary. 
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Chapter  1 

SOIL  MOISTURE  PROGRAM 


1.1  Automated  Products 


The  Soil  Moisture  program  produces  daily  temperature  and  preciptation  analyses 
and  IC-day  (“decade")  soil  moisture  estimates.  The  decade  programs  update  the  water 
balance  as  of  the  10th,  20th,  and  last  day  of  each  month.  The  third  and  final  decade 
of  each  month  also  produces  monthly  summaries. 

The  Soil  Moisture  program  covers  two  geographical  areas.  Area  1,  the  "European" 
or  "Soviet"  area,  is  bounded  roughly  by  45°N-60°fI  and  10°E-100°E.  Area  2,  the 
"Asian"  or  "China"  area  runs  from  10°N  to  5Q<>N  and  from  100°E  to  the  east  coast  of 
the  continent.  Computer  analyses  are  done  on  the  1/3-mesh  grid  (grid  distance  1/3  of 
AFGWC  full  mesh,  approximately  60  NM )  on  a  1:15,000,000  polar  stereographic  map.  The 
areas  covered  are  shown  in  Figures  1-1  and  1-2. 

In  addition  to  the  grid  point  maps,  the  decade  run  produces  area  averages  in 
tabular  form.  These  tables  list  decade  and  monthly  temperatures,  prec i pi  tat i on ,  and 
departures  from  long-term  means  for  35  regions  in  Area  1  and  for  23  regions  in  Area 
2.  For  Area  1,  the  "simple  aggregate  means"  are  simply  the  total  of  the  grid  point 
values  divided  by  the  number  of  grid  points  in  the  region.  In  Area  2,  due  to  local 
variations  of  soil,  terrain,  or  climate,  there  are  often  sparsely  sown  and  heavily 
cultivated  areas  within  a  region.  To  better  reflect  conditions  over  sown  acreage, 
grid-point  values  are  weighted  according  to  agricultural  land  use  to  produce  weighted 
means  and  departures  in  Area  2.  A  complete  list  of  the  Soil  Moisture  programs'  naps 
and  charts  foil ows : 


Daily 

Maps. 

(1) 

Maximum 

Temperature  (°C) 

(2) 

Minimum 

Temperature  (°C) 

(3) 

24-Hour 

Precipitation  (mm) 

Decade  Maps. 

(1) 

Days  in 

Decade  when  Mean  Values  Used  for  Missing 

T  enperaturcs 

(2) 

Days  in 

Decade  when  Mean  Values  Used  for  Missing 

Prec i pi t  at i on 

(3) 

Mean  Temperature  (tenths  of  °C ) 

(4) 

Decade 

Total  Precipitation  (mm) 

(  b  ) 

So i 1  Moi st ure  (mm ) 

(6) 

Ratio: 

Actual  t va pot ra nsp i ra t i on/Pot ent i a  1  E vapot ra nspi ra t i on 

(7) 

Moisture  Surplus  (mm) 

(8) 

Tract i onabi 1 i ty  ( Thornthwa i t c  Indices) 

(9) 

Decade 

Departure  from  Mean  Temperature  (tenths  of 

°C  ) 

(  10  ) 

Decade 

Percent  of  Mean  Prec i p i t a t i on 

(ID 

Percent 

of  Mean  Soil  Moisture 

2 


(1)  Monthly  Total  Precipitation  (no) 

(2)  Monthly  Mean  Temperature  (tenths  of  °C ) 

(3)  Monthly  Percent  of  Mean  Precipitation 

(4)  Monthly  Departure  from  Mean  Temperature  (tenths  of  °C) 

(5)  Ratio:  Monthly  Actual  E vapot ranspi rat i on/ Pot ent i a  1  Evapotranspi ration 

(6)  Accumulated  Monthly  Total  Precipitation  since  last  April  (mm) 

(7)  Accumulated  Monthly  Total  Precipitation  since  last  September  (nm) 

(8)  Percent  of  Monthly  Accumulated  Precipitation  beginning  in  April 

(9)  Percent  of  Normal  Monthly  Accumulated  Precipitation  beginning  in 
September 

d •  Decade  Area  Average  Tables. 

(1)  Mean  Temperature  (tenths  of  0 C ) 

(2)  Total  Precipitation  (tenths  of  mm) 

(3)  Soil  Moisture  (tenths  of  mm) 

(4)  Decade  Departure  from  Long-Term  Mean  Temperature  (tenths  of  °C) 

(5)  Percent  of  Mean  Precipitation 

(6)  Percent  of  Mean  Soil  Moisture 

e.  Monthly  Area  Average  Tables. 

(1)  Total  Monthly  Precipitation  (tenths  of  mm) 

(2)  Monthly  Mean  Temperature  (tenths  of  °C) 

(3)  Departure  from  Mean  Monthly  Temperature  (tenths  of  °C) 

(4)  Percent  of  Mean  Monthly  Precipitation 

(5)  Total  Precipitation  beginning  last  April  (tenths  of  mm) 

(6)  Total  Precipitation  beginning  last  September  (tenths  of  nr.i) 

(7)  Percent  of  Normal  Accumulated  Precipitation  since  last  April 

(8)  Percent  of  Normal  Accumulated  Preci pi  tat i on  since  last  September 

The  daily  temperature  and  precipitation  analyses  are  normally  completed  within 
two  days  after  the  data  date.  The  decade  and  monthly  products  are  normally  completed 
within  two  days  after  the  last  day  of  the  decade  period. 

1 .2  Water  Balance.  Bookkeeping  Method 

1.2.1  Water  Balance  Equation.  Soil  moisture  is  calculated  as  "the  balance  between 
the  income  of  water  by  precipitation  and  the  outflow  of  water  by  evapotranspi rati  on 
(Thornthwalte  and  Mather,  1957).  The  bookkeeping  method  consists  of  periodically 
adding  rainfall  and  subtracting  evapotranspi rat i on  and  water  surplus  to  compute  a  new 
soil  moisture  balance.  The  water  balance  equation  may  be  written 


where 


SM  =  new  soil  moisture 

SM0  =  starting  soil  moisture 

R  =  ral nf al  1 

El  »  evapot ransp i rat i on 

PERC  =  percolation 

RO  =  runoff 

According  to  Thornthwai te ,  water  in  the  process  of  percolation  is  detained  in  the 
soil  only  briefly.  "For  instance,  in  making  monthly  rather  than  daily  computations 
it  is  possible  to  consider  surface  runoff  and  the  percolation  of  gravitational  water 
as  one  quantity,  moisture  surplus"  (Thornthwaite  and  Mather,  1955).  The  water 
balance  equation  then  reduces  to 

SM  =  SM0  +  R  -  ET  -  S  (mm)  (2) 

where  S  =  PERC  +  RO  =  water  surplus. 

On  the  right-hand  side  of  Equation  (2),  SH0  is  known  from  previous  analysis,  R  is 
obtained  from  surface  observations,  ET  is  estimated  from  potential  evapotranspi rati  on 
using  Thornthwaite 's  equation,  and  S  is  any  amount  in  excess  of  field  capacity.  5M0, 
R,  and  S  are  relatively  easy  to  obtain.  Estimating  ET,  however,  first  requires  an 
estimate  of  the  potential  evapotranspi rati  on  and  is  the  hydrometeorological  heart  of 
Thornthwaite 's  method. 

1.2.2  Thornth waite's  Evapotranspiration  Formula.  Thornthwa i te  1  s  empirically  derived 
equation  for  potential  evapotranspiration  is 

PET  =  DC  1  -  6  (~i)A]  (3) 

where  D  =  mean  duration  of  daylight  (hrs) 

T  =  mean  daily  temperature  for  the  calculation  period  (°C) 

1 2 

I  =  j  (Tk/5)1'^4  is  Thornthwaites  heat  index  ( T ^  =  mean  monthly 
temperature  ) 

A  =  0.0000006751  3  -  0.0000771I2  «■  0.017921  +  0.49239 

The  value  of  I  for  a  place  is  established  from  climatology,  and  values  of  A  have 
been  tabulated  as  a  function  of  1  for  easy  application  of  the  equation.  Potential 
evapotranspiration  (PET)  in  Equation  (3)  is  also  called  "actual  PET"  by  Thornthwaite 
to  distinguish  from  "unadjusted  PET"  (PET  before  the  daylight  factor  D  is  applied) 
which  is  the  quantity  in  brackets  in  Equation  (3). 

The  handy  thing  about  Tho rnt hwa i t e 1  s  equation  is  that,  once  A,  1,  and  D  have  been 
established  for  a  site,  the  only  meteorological  element  needed  is  mean  daily 
temperature . 

In  refining  this  formula  to  achieve  more  satisfactory  results  and 
to  utilize  only  climatic  data  which  are  generally  available  it  was 
possible  to  eliminate  all  factors  but  mean  temperature  and  average 
length  of  day.  1 nat  satisfactory  results  could  be  obtained  without 
the  use  of  wind,  humidity,  or  solar  radiation  seems  to  be  due  to  the 
fact  that  all  these  important  influences  on  evaporation  including 
temperature  vary  together  (Jhorthwaite  and  Mather,  1955). 

Thornthwaite  also  found  that  mean  daily  temperature  was  the  most  suitable 
meteorological  variable  to  use  for  theoretical  reasons. 

In  developing  a  formula  which  uses  climatic  factors  for  computing 
potential  evapotranspiration  a  conservative  climatic  parameter  must 
be  used;  one  that  will  be  relatively  unaffected  by  the  introduction 
of  conditions  of  potential  evapot ra nsp i ra t i on.  Atmospheric  moisture 
is  very  sensitive  to  an  increase  of  soil  moisture  necessary  for 
potential  evapot ranspi rat  i  on  and  is  thus  unsuitable.  Maximum  and 
minimum  temperatures  are  also  both  affected;  maximum  temperature  is 
not  as  high  over  moist  soil  and  minimum  temperature  is  not  as  low. 


Thus  diurnal  range  of  temperature  is  not  a  conservative  property 
either  but  would  exhibit  a  reduction  if  the  soil  became  moist.  Since 
maximum  and  minimum  daily  temperatures  are  affected  in  the  opposite 
direction  by  changes  in  soil  moisture,  the  mean  is  only  slightly 
affected;  mean  daily  temperature  is  one  of  the  most  conservative 
climatic  elements.  Temperature  can  serve  as  an  index  to  potential 
e vapot ranspi rat i on  because  there  is  a  fixed  relation  between  the  net 
radiation  used  for  heating  and  that  used  for  evaporation  when 
conditions  exist  to  achieve  the  potential  rate  (Thornthwai te  and 
Mather,  1955). 

Once  PET  has  been  estimated,  an  estimate  of  ET  is  derived  from  it  for  use  in  the 
water  balance  equation.  If  the  soil  moisture  is  at  field  capacity,  ET  =  PtT. 
However,  as  the  soil  dries,  additional  evapot ransp i rat i on  takes  place  at  an 
increasingly  lower  rate.  Thornthwaite ‘s  studies  indicated  that  the  rate  of 
evapot rans pi  rat i on  is  proportional  to  the  amount  of  water  remaining  in  the  soil.  For 
example,  if  the  soil  moisture  is  at  one-quarter  of  field  capacity,  ET  =  PET/4. 

Thorthwai te '  s  method  has  been  applied  in  varied  regions  with  varied  results.  In 
general,  it  works  best  in  mid-latitude  climates  in  which  it  was  developed  and  gives 
better  results  when  applied  over  longer  time  periods  (Chang,  1968). 

1 . 3  Soil  Moisture  Production  Programs 

1.3.1  Daily  Analysis.  The  surface  data  for  the  Soil  Moisture  temperature  and 
precipitation  analyses  are  collected  at  AFCWC  via  the  Automated  Weather  Network  (AWN) 
and  sent  to  USAFETAC  each  day.  At  USAFL  TAC ,  the  observations  are  put  into  a 
station-order  sort  and  saved  as  part  of  USAFETAL's  climatological  data  base.  The 
Soil  Moisture  analysis  nrograms  are  run  immediately  after  the  daily  surface  sore 
program.  Over  15,000  3-hourly  surface  synoptic  observations  for  the  Soil  Moisture 
areas  are  received  each  day.  Of  these,  some  5,500  observations  contain 
maxi  mum/mi ni nun  temperature  and/or  precipitation  reports. 

The  Soil  Moisture  daily  analysis  is  done  in  two  job  steps.  The  first  step 
extracts  temperature  and  precipitation  data,  assigns  the  data  to  grid  points,  and 
removes  duplicate  reports.  The  second  step  analyzes  the  data  over  the  two  Soil 
Moisture  areas.  Figure  1-3  shows  the  daily  Soil  Moisture  production  stream. 


Figure  1-3.  Daily  Soil  Moisture  Run  Stream. 


1.3.2  Analysis  Method.  The  Soil  Moisture  temperature  and  precipitation  analyses  are 
done  using  the  least-squares  surface  fitting  technique  of  Mount,  et  al.  (1962).  The 
object  of  the  least-squares  method  is  to  minimize  the  difference  between  analysis 
values  and  actual  values.  The  Soil  Moisture  program  uses  the  1  east -squares  nethod  to 
fit  a  plane  surface  to  the  observations  around  each  grid  point;  the  low  density  of 
observations  over  most  of  the  Soil  Moisture  grid  prohibits  use  of  a  more  conplex 
quadratic  surface.  The  observed-analysis  difference  to  be  minimized  is: 

e  (Ho  -  Ha)?  (4) 

where  H0  =  observed  value 
Ha  =  analysis  value 

A  plane  through  the  point  Ha  is  given  by 

Ha  =  Ax  +  By  +  C  (S) 


As  each  grid  point  is  analyzed,  it's  coordinates  are  taken  as  (0,0)  so  that  C 
determines  the  analysis  value,  Ha  (0,0)  =  C.  The  x  and  y  coordinates  are  then  taken 
relative  to  the  analysis  point. 

e  is  minimized  when 

ff  -  z  [x  (H0  -  Ha)]  =  0, 

If  =  z  C  y  (H0  -  Ha ) ]  =  0,  and 

|f=l(Ho-Ha)=0  (6) 

Substituting  Equation  (5)  into  Equation  (6)  and  rearranging  gives 
ZH0x  =  zAx2  +  ZB xy  +  ZCx, 

Zll0y  -  zAxy  +  ZBy2  +  zCy,  and 

ZHQ  =  eAx  +  ZBy  +  C  •  N  (  7 

where  N  is  the  number  of  observations. 

Applying  Cramer's  Rule  (see,  for  example,  Korn  and  Korn,  1^61)  to  Equation  (7),  we 
get 


C  -  Ha  (0,0)  = 


zx2 

Zxy 

X 

O 

X 

txy 

zy  2 

zH0y 

LX 

zy 

zh0  i 

L  X  ^ 

Zxy 

ZX  j 

zxy 

CM 

>> 

13 

zy  | 

LX 

zy 

N  ! 

;H0y: 

:xy-: 

:li0xzy2) 

zx(zyzxy-zxzy2)+zx2(Nzy2-zyzy)+zxy(zxzy-Nzxy) 
where  C  is  ehe  analysis  value. 
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At  least  three  observations  are  needed  to  uniquely  determine  a  plane  through  the 
data.  On  its  first  data  scan,  the  Soil  Moisture  program  considers  only  observations 
within  one  grid  distance  of  the  analysis  point.  If  there  are  three  or  more 
observations,  the  least-squares  plane  value  is  used.  If  only  two  observations  are 
available,  a  simple  mean  value  is  used.  If  there  is  only  one  observation  and  it  is 
within  1/2-grid  distance,  that  observation  value  is  used.  If  a  value  cannot  be 
assigned,  the  program  goes  out  two  grid  distances  for  the  precipitation  analysis  and 
up  to  four  grid  distances  for  the  temperature  analyses,  searching  for  three  or  more 
observations.  If  a  plane  still  cannot  be  determined,  the  analysis  value  is  filled 
with  a  missing  indicator. 

The  daily  analyses  are  quality  controlled  by  USAFETAC  analysts.  Missing  and 
suspect  analysis  points  are  checked  using  data  printouts,  nearby  reports,  and  spatial 
continuity.  Corrections  to  the  analyses  are  made  by  card  input  to  the  Soi  1  Moisture 
bookkeeping  programs.  If  a  value  is  left  missing,  the  bookkeeping  programs  insert  a 
mean  climatological  value. 

1.3.3  Soil  Moisture  Decade  Production.  The  soil  moisture  bookkeeping  parameters  are 
updated  three  times  a  month,  following  the  daily  processing  for  the  10th,  20th,  and 
last  day  of  each  month,  in  what  is  known  as  the  decade  run.  The  decade  prograns  use 
the  daily  precipitation  and  temperature  analyses  to  estinate  the  terms  needed  to 
solve  the  water  balance  equation  as  described  in  Section  1.2.1.  The  calculated  water 
balance  quantities  are  used  to  update  the  soil  moisture  bookkeeping  from  the  last 
decade.  The  decade  run  consists  of  three  job  steps  as  shown  in  Figure  1-4. 
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Figure  1-4.  Soil  Moisture  Decade  Run. 
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Figure  1-5.  Soil  Moisture  Calculation. 
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The  first  decade  job  merges  the  dally  Soil  Moisture  tapes  and  Inserts  analysis 
corrections  fron  card  Input. 

The  second  step  does  the  water  balance  calculations,  updates  the  bookkeeping 
parameters,  and  prints  the  grldded  maps.  Figure  1-5  shows  the  major  steps  in  the 

soil  moisture  calculations.  Note  the  Soil  Moisture  assumption  of  200-mn  field 
capacity  at  all  points.  Soil  moisture  values  over  200  are  set  to  200  and  the 

remainder  Is  counted  as  water  surplus. 

The  third  decade  job  prepares  the  district  averages  In  tabular  form  (these  a  re 
described  in  Section  1.1)  from  the  updated  bookkeeping  parameters. 

1 . 4  Summary 

The  Soil  Moisture  program  is  limited  by  its  many  simplifying  assumptions.  These 
limitations  include 

a.  The  Thornthwaite  evapot ranspi rat i on  formula  does  not  consider  turbulent 

transfer  or  thermal  advection.  These  are  important  effects  on  ET,  especially  on 

shorter  time  scales  (Thornthwaite  and  Hare,  1965;  Chang,  1968). 

b.  The  assumption  of  200-nn  field  capacity  at  all  points. 

c.  ET  is  set  to  zero  if  the  mean  decade  temperature  is  less  than  0°C. 

d.  The  data  are  sparse  over  some  parts  of  the  grid. 

e.  The  bookkeeping  inherently  assumes  that  precipitation  is  spread  evenly  over 
the  decade  period.  If  most  of  the  precipitation  actually  occurs  early  in  the  decade, 
actual  soil  moisture  will  be  less  than  the  calculated  value,  and  vice  versa. 

These  simplifications,  of  course,  make  the  Soil  Moisture  program  easy  to  use. 
Based  on  Thornthwaite  and  Mather's  (1957)  instructions  for  soil  moisture 
calculations,  a  major  objective  of  the  method  was  that  it  required  only  routinely 
available  surface  data  as  input.  These  observations  do  not  give  the  information 
needed  for  a  more  detailed  treatment.  However,  the  Thornthwaite  method  has  been 
found  to  give  good  results  in  many  applications.  The  U5AFETAC  Soil  Moisture  program, 
with  it's  10-day  time  intervals,  probably  uses  as  small  a  time  resolution  as  is 
feasible.  The  Soil  Moisture  program  is  simple  and  easy  to  use,  and  it  has  provided 
useful  results  for  over  20  years. 


Chapter ; 2 
AGROMET  PROGRAM 


2.1  Automated  Products 


The  Agromet  program  produces  daily  grid-point  values  of  evaporation  from  a  free 
water  surface  (called  evapotranspi ration  potential  or  ETP  in  the  program),  the 
meteorological  elements  used  to  estimate  ETP,  and  a  modeled  stage  of  plant  growth 
(phenology).  The  phenology  and  analyses  of  maximum  temperature,  minimum  temperature, 
precipitation,  and  snow  depth  are  produced  in  nap  form.  ETP,  net  solar  radiation, 
and  accumulated  daily  wind  are  listed  by  grid  point.  Agronet  uses  the  1/8-ncsh  AFGWC 
grid  on  a  polar  stereographic  projection,  which  gives  approximately  25-NM  grid  point 
spacing  at  60°N. 

The  Agromet  program  covers  six  geographical  areas.  Area  1  includes  eastern 
Europe  and  the  central  USSR.  Area  2  covers  China  east  of  100°E.  Area  3  covers  the 
Bay  of  Bengal  region  and  includes  Burma,  Bangladesh,  and  NE  India.  Area  4  covers  NW 
India,  Pakistan,  Afghanistan  and  Iran.  Area  5  covers  most  of  the  United  States. 
Area  6  covers  Europe  east  from  10°E  to  the  boundary  of  Area  1.  The  six  Agromet  areas 
are  shown  in  Figure  2-1. 

2.2  Methods 


2.2.1  The  Penman  Equation.  Penman  (1956)  combined  estimates  of  the  radiation  and 
turbulent  diffusion  effects  on  evaporation1  in  a  single  equation.  This  was  a 
significant  improvement  over  earlier  formulations,  such  as  Thornthwai te 1 s  (Section 
1.2.2),  that  modeled  energy  requirements  only.  In  his  development,  Penman  assumed 

equality  of  the  coefficients  of  vapor  transfer  equations  (Penman,  1956;  Thornthwaite 

and  Hare,  1965).  This  assumption  allows  the  transport  constants  to  fall  out  of  the 
Penman  equation 

ftp  =  A  rN  —  ,.Y  r — IDE —  a  (g) 

t a  +  Y  L  day  J  m 

where  ETP  =  evaporation  from  an  open  water  surface 


R u  =  net  radiation  expressed  in  evaporation  units 

Ea  =  aerodynamic  component  expressed  in  evaporation  units 

A  =  d(cs)/dT,  the  slope  of  the  saturation  vapor  pressure  vs  temperature 

curve 


=  the  psychometric  constant,  ratio  of  the  specific  heat  of  dry  air  to 
the  latent  heat  of  vaporization  of  water 


The  relative  importance  of  the  radiation  and  the  aerodynamic  effects  in  Equation 
(9)  is  given  by  the  ratio  a/y.  In  the  Agroment  program,  Y  is  taken  as  constant  at 
0.64  (mb/K),  although  it  does  vary  slightly  with  temperature  and  pressure.  A  is 
calculated  from  the  equation 


A 


dc_s  « 
dT 


(6790. 5-5. 02808T+49 16. 8x10-- 03C4TT2+1 74209x10- 1302. B8/T ) 


(10) 


given  on  page  37  2,  Smithsonian  Meteorological  Tables  (List, 1951).  a  ranges  from  0,44 
at  0°C  to  3.9  at  40°C.  At  7 0  C ,  A  ■  0.64  »  Y,  and  at  18°C,  a  ■  2Y.  Thus,  at 
typical  growing  season  temperatures,  A  is  larger  than  Y  .  In  addition,  Rn  is  usually 
larger  than  Ea,  so  that  the  radiation  term  is  usually  most  important  in  Equation  (9). 
However,  the  aerodynamic  term  shows  large  short  period  and  local  variation 
(Thornthwaite  and  Hare,  1965).  By  including  the  aerodynamic  component,  Penman's 
equation  provides  more  detail  than  earlier  methods  and  is  more  complete  physically. 


Fiijirp  <>-la.  Arjro"iet  Areas  1  and  6. 


Figure  2-ld.  Aqromet  Area  4. 


Figure  2-le.  Agromet  Area 


2.2.2  The  Aerodynamic  Term.  The  aerodynamic  tern  in  Penman's  tquation  (9)  Is 
calculated  usi ng  Penman  ’ s  empirical  formula 

Ea  *  0.35  (0.5  +  yjfo)  (es  -  e)  C-ggy]  (11) 

where  u  =  cumulative  surface  wind  over  a  24-hour  period  (miles/day) 
es  »  saturation  vapor  pressure  (mb) 
e  *  vapor  pressure  (mb) 

Since  e  and  es  are  related  to  temperature,  Ea  can  be  calculated  from  temperature 
and  wind  analyses.  The  saturation  vapor  pressure  is  calculated  from  the  Goff-Gratch 
formula  given  on  page  350,  Smithsonian  Meteorological  Tables  (Listj 1951 ) 

log  io  es  =  -7.90298  £(TS/T)  -  1]  +  5.02808  1 og i0  ( Ts/T ) 

-1.3816  x  IO"9  (10H-334[1  -  (T/Ts)]  .  i) 

+8.1328  x  (io-3.49149[(Ts/T)  -  1]  .  i) 

+  logio  ews  (12) 

where  Ts  =  steam  point  temperature  =  3 7 3 ( K ) 

T  ■  air  temperature  (K) 

ews  =  saturation  vapor  pressure  at  Ts  =  1013(mb) 

The  vapor  pressure,  e,  also  needed  in  Equation  (11),  is  calculated  from  Equation 
(12)  using  dew-point  temperature,  T d ,  in  place  of  T. 

2.2.3  The  Radiation  Term.  The  net  radiation,  Rn,  needed  in  the  Penman  Equation  (9) 
is  estimated  by  considering  the  radiation  balance 

RN  =  RSE  -  Rst  +  RNL  [-^-n]  (13) 

where  Rs+  =  incoming  solar  radiation 

Rs+  =  solar  radiation  reflected  at  the  surface 

R N l,  =  net  long-wave  radiation 

The  three  terms  on  the  right-hand  side  of  Equation  (13)  are  estimated  using 
incoming  solar  radiation  calculations  and  detailed  cloud  analyses.  The  cloud 
analyses  are  used  to  model  the  incoming  solar  radiation  and  net  infrared  radiation  in 
three-hour  steps.  Albedo  is  then  estimated  as  a  function  of  crop  stage  of  growth  and 
used  to  calculate  the  reflected  solar  radiation.  These  calculations  are  described  in 
detail  below. 

a.  Direct  Beam  Solar  Radi ati on.  The  incoming  flux  of  direct  solar  radiation 
under  clear  skies  is  calculated  from  a  well-known  set  of  equations  (see,  for  example, 
Haltiner  and  Martin,  1957  ;  List, 1951).  The  basic  equation  is 

roirect  ■  J0  *ssc  !  '«  *  1141 

where  J0  «  solar  constant  =  1.94  (cal/cm2min) 

a  -  atmospheric  transmissions  coefficient  for  solar  radiation 
z  *  zenith  distance  of  the  sun 
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The  zenith  distance,  defined  as  the  angular  distance  of  the  sun  from  directly 
overhead,  is  given  by 

cos  z  =  sin  jisin  i  +cos  (  cos  '  cosh 
where  't>  =  latitude 

f,  =  solar  declination  (latitude  at  which  the  sun  is  overhead) 

h  «  hour  angle  (angular  distance  of  rotation  from  solar  noon } 

The  transmission  coefficient  for  solar  radiation  is  given  by 

a  =  0.73  -  0.077  u0-3  (lb) 

where  u  is  the  precipitable  water  vapor  (cm). 

The  factor  0.73  is  used  to  account  for  depletion  of  solar  radiation  by  dry  air, 
including  dust.  The  term  0.077  u^*3  Was  developed  by  MacDonald  (  1960)  to  account  f o  •• 
attenuation  of  solar  radiation  by  water  vapor.  The  atmospheric  water  vapor  is 

estimated  from  the  surface  vapor  pressure  using  an  empirical  relationship  developed 

by  Idso  (1969) 


logjQ  u  =  -U.  579  +  0. 274  /e  (  1 6  ' 

where  e  is  the  surface  vapor  pressure. 

From  Equations  (14),  (15),  and  (16),  surface  vapor  pressure,  latitude,  time 

year  (for  solar  declination),  and  time  of  day  (for  hour  angle)  are  the  basic  cata 
needed  to  estimate  direct  beam  solar  radiation  under  clear  skies. 


b.  Oiffuse  Solar  Radiation.  Following  the  procedure  given  in  the  Pmitnsoniai; 
Meteorol ogi cal  Tables  (List,  T951,  p.  420),  diffuse  solar  radiation  is  calculated  as 


RD T  FFUSE 


0.91  R  - 
_ o 

2 


DIRECT 


cm 


mi  n 


where  R0  =  cos  z  is  the  extra-terrestrial  incoming  solar  radiation  flux.  Iron 

Equations  (14)  and  (17),  the  total  solar  radiation  reaching  the  surface  under  clear 
skies  is 

4-  +  + 

rCUAR  =  rDIKECT  +  rD  I FFUPE  (ic 


c.  Polar  Radiation  Cloud  Factors.  The  Agromet  program  uses  the  AFGWC  Automated 
Cloud  Analysis  Model  ( 3DNEPH )  data  at  1/8-mesh  grid  points  (see  Fye,  19/B)  to  adjust 
the  clear  sky  solar  radiation  for  cloudiness.  Table  2-1  shows  the  15  cloud  layers 
used  in  the  30NEPH  model.  30NEPH  provides  percent  cloud  cover  in  each  layer,  total 
cloud  amount,  and  low,  middle,  and  high  cloud  types  (Table  2-2). 
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Table  2-1.  3DNEPH  Cloud  Layers. 


LAYER  NO. 

LAYER 

INTERVAL 

1 

SFC  - 

150  FT  AGL 

2 

151  - 

300  FT  AGL 

3 

301  - 

600  FT  AGL 

4 

601  - 

1000  FT  AGL 

5 

1001  - 

2000  FT  AGL 

6 

2001  - 

3500  FT  AGL 

7 

350’.  - 

5000  FT  MSL 

8 

5001  - 

6500  FT  MSL 

9 

6501  - 

10000  FT  MSL 

iO 

10001  - 

14000  FT  MSL 

1  1 

14001  - 

18000  FT  MSL 

12 

18001  - 

22000  FT  MSL 

13 

22001  - 

26000  FT  MSL 

14 

26001  - 

35000  FT  MSL 

15 

35001  - 

ABOVE  FT  MSL 

Table  2-2.  30NEPH  Cloud  Types. 

CODE 

LOW  CLOUD  TYPE  3 

MIDDLE  CLOUD  TYPES 

0 

Type  Unknown 

Type  Unknown 

1 

Stratocunulus  (SC) 

Altocumulus  (AC) 

2 

Stratus  (ST) 

Altostratus  (AS) 

3 

Curnu  1  u s  ( CO  ) 

Ninbostratus  (NS) 

4 

Cumulonimbus  (CB) 

AC  and  AS 

5 

SC  and  ST 

AC  and  NS 

6 

SC  and  CU 

AS  and  NS 

7 

SC  and  CB 

AC  and  AS  and  NS 

8 

ST  and  CU 

9 

ST  and  CB 

10 

CU  and  CB 

11 

SC  and  ST  and  CU 

12 

SC  and  ST  and  CB 

13 

SC  and  CU  and  CB 

14 

ST  and  CU  and  CB 

1  5 

SC  and  ST  and  i U  and 

i  B 
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HIGH  CLOUD  TYPE 
Type  Unknown 
Cirrus  (Cl) 
Cirrocunulus  i  •  ■ 
Cirrostratus  ( c  S  , 
Cl  and  C C 
Cl  and  C  S 
C  C  and  C  S 
cl  and  C C  and  LG 


■fiBBBBnwMdJ 


Agromet  reduces  3DNEPH  cloud  cover  to  percentage  cloud  cover  in  each  of  three 
layers:  low,  middle,  and  high.  A  "solar  radiation  cloud  coefficient"  (see  Smithso¬ 
nian  Meteorological  Tables,  p.  441  of  List,  1951),  is  associated  with  each  type  in 
each  of  the  three  layers  (Table  2-3). 


Table  2-3.  Solar 

Radiation 

Cloud 

Coefficient  (Ks) 

as  a  Function  of  Cloud 

Solar 

Radiation  Cloud 

Coefficient  (Ks) 

30NEPH 

Cloud  Type* 

Low 

Cloud 

Middle  Cloud  High  Cloud 

0 

1.0 

1.0 

1.0 

1 

0.34 

0.51 

0.82 

2 

0.25 

0.41 

0.82 

3 

0.34 

0.17 

0.82 

4 

0.12 

0.46 

0.82 

5 

0.30 

0.30 

0.82 

6 

0.34 

0.25 

0.82 

7 

0.20 

0.35 

0.82 

8 

0.30 

1.0 

1.0 

9 

0.18 

1.0 

1.0 

10 

0.20 

1.0 

1.0 

11 

0.32 

1.0 

1.0 

12 

0.22 

1.0 

1.0 

13 

0.24 

1.0 

1.0 

14 

0.22 

1.0 

1.0 

15 

0.26 

1.0 

1.0 

♦See  Table  2-2. 


Using  these  cloud  coefficients  and  the  low,  middle,  high  cloud  amounts  calculated 
from  the  3DNEPH,  a  cloud  solar  radiation  depletion  factor  is  calculated 


3 

F  =  n  [1  -  C .  ( 1  -  K  )] 
s  i  =  l  1  s  i 


(19) 


where 

Cl 

*  the  cloud  amount 

1 

=  1,  2,  3  for  low,  middle,  and  high  cloud 

KS 

*  the  solar  radiation  cloud  coefficient  for  the 
3DNEPH  In  layer  i 

cloud  type  reported  by 

Fs 

■  the  fraction  of  the  Incoming  solar  radiation 
cover 

depleted  due  to  cloud 
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Combining  Equations  (18)  and  (19),  the  Incoming  solar  radiation  Is  given  by 


Fs  rsclEAR 


r  C  4>  1  *i 

cm^ml n 


(20) 


by 


d.  Reflected  Solar  Radiation.  Solar  radiation  reflected  at  the  surface  Is  given 


Rs ♦  «  -A  Rs+ 


r  Cd  1  ^ 

cm2ni  n 


(2D 


where  A  Is  the  albedo. 

In  the  Agromet  program,  the  albedo  Is  specified  as  a  function  of  the  modeled 
stage  of  crop  growth  (phenology).  The  relationship  is  shown  In  Table  2-4. 


Table  2-4.  Richardson's  Relationshop 
Phenology 
-3.0  <  P  <  0.99 
1.0  <  P  <  1.99 
2.0  <  P  <  3.99 
4.0  <  P  <  5.99 
6.0  <  P  <  6.99 


for  Albedo  as  a  Function  of  Phenology. 

0.90 

0.90  -  0.14  (P  -  1.0) 

0.76 

0.76  +  0.14  (P  -  4.0) 

0.90 


Combining  Equations  (20)  and  (21),  the  net  solar  radiation  is 

RNS  -  Rs+  +  M  *  (1-A)  RST  Cc(4^n . 3  (22) 

e.  Long-Have  Radiation.  The  net  long-wave  radiation  is  calculated  following  Geiger  (1965). 

For  clear  skies 

RN1  -  cr  T4  (0.18+0. 25X10-°*126e)-0.007(T-Ta)  [— §-■  -]  (23) 

CLEAR  cmsni  n 

where  a  =  Stephan-Boltzmann  constant 
T  =  air  temperature  (K) 

Tg  =  ground  temperature  (K) 
e  *  vapor  pressure  (na  Hg) 

In  the  Agromet  program,  T  «  Tg  Is  assumed,  and  e  is  in  mb.  Equation  (23)  reduces  to 


r  cal  *1 
cra^ml n 


RNLcLEAR  “  a  t4<0-18  +  0.25  *  io-0.0945ej 


(24) 


The  effects  of  clouds  are  accounted  for  by  a  long-wave  radiation  cloud  factor,  Fl 
(Geiger,  1965).  The  method  used  to  get  F|_  is  similar  to  that  used  for  the  cloud 
solar  radiation  depletion  factor  (Section  2.2.3c).  Individual  long-wave  cloud 
coefficients  ( Kl )  are  associated  with  the  cloud  type  in  each  of  the  low,  middle,  and 
high  cloud  layers  (Table  2-5).  Fl  is  calculated  by 


FL  =  {  2  [C./iTT  *  (1  -  C.)]}: 

i=l  1  L1  j=o  J 


where 


cloud  amount  in  layer  i  (C0  =  0  by  definition) 
i^l  for  the  cloud  type  in  layer  i 


Table  2-5.  Long-Uave  Cloud  Coefficient  (Kl)  as  a  Function  of  Cloud  Type  (3DNEPH). 


3DNEPH  Code 


Long-Wave  Coefficient  (Kl) 


All  Types 

£l 

Cm 

C_H 

0 

0 

0 

0 

1 

0.20 

0.17 

0.04 

2 

0.24 

0.20 

0.06 

3 

0.20 

0.24 

0.08 

4 

0.24 

0.19 

0.06 

5 

0.22 

0.21 

0.06 

6 

0.20 

0.21 

0.06 

7 

0.22 

0.21 

0.06 

8 

0.22 

0 

0 

9 

0.24 

0 

0 

10 

0.22 

0 

0 

11 

0.22 

0 

0 

12 

0.22 

0 

0 

13 

0.22 

0 

0 

14 

0.22 

0 

0 

15 

0.22 

0 

0 

1 ong-wave 

radiation,  including  cloud  effects. 

is  then 

given 

rnl 

*  Rjkj»  -  F.  (o  T  -  Rml  ) 

nlclear  l  "lclear 

j-  cal 
cm2ni n 

] 

2.2.4  Agromet  ETP  Calculation.  With  the  estimates  of  the  three  radiation 
components ,  Rf ,  Rf,  and  R^l»  and  of  the  aerodynamic  term,  E»,  the  evapotranspi rat i on 
potential  nay  be  calculated  from  the  Penman  Equation  (9).  In  practice,  the 
evapotranspi rat i on  from  the  incoming  solar  radiation,  ETP1,  and  that  from  long-wave 
radiation  and  the  aerodynamic  term,  ETP2,  are  calaculated  seperately  at  AFGWC  for 
each  grid  point.  At  USAFETAC,  the  temperature  and  precipitation  analyses  are  used  to 
run  the  plant  phenology  model,  and  albedo  is  estimated  as  a  function  of  plant 
phenology.  ETP  is  then  calculated  from  ETP1,  ETP2,  and  albedo. 

The  form  of  the  Penman  equation  used  by  Ayromet  is  developed  from  substituting 
Equation  (13)  into  Equation  (9) 


ETP  = 


A  (Rst  4  RsQ  +  A  Rn|_  ♦  y  Ea 
5  +  y  A  +  y 


(27) 


Then,  substituting  from  Equation  (21) 

i_V  ' 

A  +  y  A  +  y 

And,  by  definition  of  ETP1  +  ETP2 


s'  ’  R,.I  +  Y  E 

ETP  =  (1  -  A)  — +  —Jit - a 


(23) 


ETP  =  (1  -  A)  ETP1  +  ETP2 


(29) 


2.3  Agromet  Production  System 

2.3.1  AFGWC  Ag  rome  t .  The  AFGWC  Agromet  production  system  runs  eight  times  a  day, 
following  the  processing  of  3-hourly  (00z  plus  every  3  hours)  3DNEPH  cloud  analysis 
model  (Fye,  1978)  and  surface  observations.  The  00Z  through  18Z  Agromet  runs  consist 
of  three  job  steps:  a  surface  analysis,  a  first  guess  for  daily  maximum  and  minimum 
temperature,  and  a  summation  of  the  Penman  equation  parameters  and  estimated  precipi¬ 
tation.  In  addition  to  these  three  steps,  the  21Z  run  produces  the  final  daily  tem¬ 
perature  and  precipitation  analyses  and  the  daily  Agromet  data  tape.  Figure  2-2 
shows  the  AFGWC  production.  The  individual  programs  are  discussed  briefly  in  the 
fol 1 ow  i  n  g  sections. 

a.  Surface  Analysis.  The  first  program  to  run  in  eacn  3-hourly  Ayromet  cycle 
is  the  surface  analysis.  This  program  provides  a  1/4-mesh  grid  point  analysis  of 
surface  temperature,  dew-point  depression,  and  u  and  v  wind  components. 

These  analyses  are  done  over  an  expanded  area  which  extends  one  full  mesh  grid 
distance  beyond  the  boundaries  of  each  Agromet  area.  This  procedure  excludes 
possible  boundary  problems  from  the  stored  analysis  fields. 

The  surface  temperature  and  dew-point  depression  analyses  use  the  Barnes  (1964) 
technique  with  three  scans.  The  u  and  v  wind  components  are  derived  from  1000-ub 
D-values.  The  D-values,  in  turn,  are  derived  from  a  Barnes  analysis  of  surface 
temperature  and  pressure.  The  latest  available  coarse  mesh  analysis  (dew  point)  and 
forecast  (temperature  and  pressure)  from  the  AFGWC  data  base  are  used  as  first  yuess 
for  these  analyses.  These  analyses  are  stored  in  the  1/4-mesh  analysis  data  base 
which  is  created  by  Agromet. 

b.  Maximum  and  Minimum  Temperature  First  Guess.  This  program  stores  daily 
maximum  and  minimum  temperatures  from  the  AFGWC  1/2-mesh  hemispheric  surface 
temperature  analysis.  Every  3  hours,  the  1/2-mesh  analysis  is  interpolated  to  1/4 
mesh  and  the  value  at  each  grid  point  compared  to  the  maximum  and  minimum  values 
found  thus  far  in  the  day.  If  the  latest  temperature  exceeds  the  previous  maximum  or 
minimum,  the  new  extreme  replaces  the  previous  value.  The  resulting  tields  are 
stored  in  the  1/4-mesh  data  base  and  are  used  as  first  yuess  by  the  21Z  cycle 
maximum/minimum  temperature  analysis. 
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Figure  2-2.  AFGWC  Daily  Agromet  Production. 
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c.  ETP  Calculations  and  Estimated  Precipitation.  This  program  calculates  the 
evapotranspi rati  on  (Etp)  parameters  ( see  Section  2T  2 )  and  estimated  precipitation 
every  3  hours.  The  tTP  parameters  are  calculated  using  the  temperature,  dew-point 
depression,  and  surface  Mind  components  stored  in  the  1/4-mesh  data  base  and  the 
AFGWC  3DNEPH  model  data  base.  Solar  radiation  reaching  the  surface  is  computed  at 
each  grid  point  by  positioning  the  sun  at  its  location  1-1/2  hours  before  data  time 
and  moving  it  across  the  sky  in  1/2-hour  increments  to  its  position  1-1/2  hours  after 
data  time.  Daily  values  of  the  ETP  parameters  thus  represent  the  period  from  2230/ 
of  the  previous  day  through  2230Z  on  the  data  date.  The  ETP  parameters  are  summed 
every  3  hours  and  stored  in  the  Agronet  data  base. 

An  estimated  prec i pi  tat i on  value  is  calculated  for  each  grid  point  every  3  hours 
buied  on  the  cloud  type  and  amount  from  the  3DNEPH  data  base  using  a  method  suggested 
by  Follansbee  (1973).  The  estimated  precipitation  is  accumulated  at  each  grid  point 
through  the  data  day. 

Three-hour  quantitative  precipitation  forecast  (OFF)  values  are  also  stored  at 
eauh  grid  point.  The  QPF  values  are  from  the  AFGWC  b-layer,  1/2-oesh  cloud  forecast 
model.  The  1/2-mesh  QPF  field  is  interpolated  to  1/8-mesh  and  the  sum  of  the  QPF 
values  at  each  grid  point  is  accumulated  during  each  Agromet  run  and  stored  in  the 
Agronet  data  base. 

d .  Max  inun  and  Minimum  Temperature  1/4-Mesh  Analysis.  The  Agronet  maximum  and 
minimum  daily  temperature  analysis  is  done  on  the  212  run.  Using  maximum  and  minimum 
temperature  reports  from  surface  observations,  the  program  does  a  Barnes  (1964) 
analysis  with  three  scans  for  both  fields.  The  first  guess  fields  are  supplied  fron 
the  3-hourly  Agromet  runs  (see  Section  2.3.1b).  These  analyses  are  stored  in  the 
1/4-mesh  analysis  data  base. 

e.  Precipitation  Total  from  Surface  Reports.  This  program  assigns  precipitation 
reports  from  surface  observations  to  the  nearest  1/8-mesh  grid  point  to  the  reporting 
station.  Two  values  are  accumulated,  precipitation  amounts  as  reported  and  "bogus" 
precipitation  values.  The  bogus  precipitation  values  are  estimates  of  precipitation 
based  on  the  weather  types  reported  in  the  observation.  These  values  are  stored  by 
the  AFGWC  data  base  building  programs  in  a  reserved  part  of  the  surface  observations. 
Agromet  uses  the  accumulated  reported  values  if  they  are  greater  than  b  mm. 
Otherwise,  it  uses  the  greater  of  the  two  values.  The  program  also  checks  for 
convective  precipitation  and  spreads  one-half  of  the  daily  accumulation  to 
immediately  adjacent  grid  points  if  no  other  precipitation  is  reported  for  those  grid 
points.  The  resulting  analysis  is  stored  in  the  Agromet  data  base. 

f.  24-Hour  Precipitation  Analysis.  This  program  spreads  the  precipitation 
tot  «• 1  •  a. cumulated  from  surface  observation  to  surrounding  grid  po'rts  with  -it. 
r-p.  .he  spread  is  based  on  the  precipitation  amounts  estimated  I'cp  the  31'Ntl’F- 
cloud  analysis.  An  estimated  precipitation  amount  is  available  at  each  grid  point. 
For  points  that  also  have  reported  precipitation  amounts,  the  ratio  Reported 
Anount/tstinated  Amount  (R/E)  is  calculated.  The  (R/E)  values  are  then  spread  using 
a  linear  distance  weighted  scheme.  Each  grid  point  then  has  an  (R/F|  value  or  a 
spread  value,  (R/E)s  associated  with  it.  Reported  prec i pi  tat  ion  is  used  tor  all 
points  where  available.  For  points  with  no  reported  precipitation,  the  (R/E. )  s  is 
used  to  determine  a  value.  If  the  (R/E )  s  is  less  than  10,  precipitation  is  (R/t)s 
times  the  estimated  value  from  the  3DNEPH;  otherwise  the  QPF  value  is  used  for  the 
grid  point.  Normally,  less  than  12  QPF  values  are  used  in  any  Agromet  area.  The 
precipitation  analysis  scheme  is  shown  in  Figure  2-3. 

g.  Agromet  Tape  Build.  The  daily  Agromet  data  tape  is  produced  on  the  21 L 
cycle.  The  program  performs  gross  error  checks  on  the  final  values  and  checks  the 
tapes  for  proper  format.  Finally,  the  program  resets  all  Agromet  data  base  fields  to 
zero  for  the  start  of  the  next  day's  run. 

2.3.2  USAFETAC  Programs.  Agromet  processing  at  USAFETAt  starts  with  the  receipt  of 
the  dally  Agromet  data  tape  from  AFGWC.  The  tape  is  checked  shortly  after  trans¬ 
mission  to  ensure  that  all  data  records  are  in  the  right  format,  that  all  characters 
are  valid,  and  that  the  temperature  and  precipitation  fields  are  filled.  The  check 
program  also  reformats  the  data  for  processing  on  USAFETAC's  IBM  4341  computer.  Unce 
a  tape  has  been  received  and  checked,  the  daily  Agromet  production  is  scheduled. 
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The  da’ly  Agronet  production  at  USAFETAC  consists  of  four  steps,  shown  in  figure 
2-4.  f'.st,  grid-point  maps  of  maximum  temperature,  minimum  temperature, 
precipitai  .on,  and  snow  depth  are  prepared  and  written  to  tape  and  hard  copy.  The 
second  step  updates  the  crop  phenology  model  based  on  the  meteorological  parameters. 
Third,  albedo  is  computed  based  on  the  phenology,  and  the  albedo  is  used  to  compute 
net  solar  radiation  and  ETP.  These  first  three  steps  are  repeated  for  each  Agronet 
area.  A  separate  file  is  written  to  the  output  tape  for  each  of  the  six  areas.  The 
fourth  step  of  the  Agronet  production  copies  observations  from  the  USAFtTAC  surface 
observation  file  onto  tne  courier  tape.  The  final  output  is  then  a  seven-file  data 
tape  for  transmission  to  Agromet  customers. 


Figure  2-3.  Precipitation  Analysis  Scheme. 
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F i  gure  2-  4 . 
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GLOSSARY  OF  SOIL  MOISTURE  TERMS 


a.  Soil  Moisture.  The  amount  of  water,  measured  as  depth  over  an  area,  retained 
against  gravity  in  the  root  zone  of  the  soil-  The  root  zone  is  the  depth  of  soil 
from  the  surface  to  the  bottom,  or  just  below  the  plant  roots.  Soil  moisture  depends 
on  moisture  availability,  soil  type,  vegetation  type,  and  plant  stage  of  growth. 

b.  Field  Capacity.  The  maximum  amount  of  soil  moisture.  Field  capacity  varies 
from  a  few  millimeters  in  shallow  sand  to  well  over  400  mm  in  a  deep  silt  loam. 

c .  Mater  Surplus.  Hater  in  excess  of  field  capacity.  This  water  cannot  be 
retained  and  is  in  the  process  of  runoff,  percolation,  or  evaporation. 

d.  Runoff .  Hater  discharged  from  an  area  through  stream  channels. 

e.  Perco 1  at i on .  The  gravi tat i ona 1  pull  of  water  from  the  root  zone  into  ground 
water.  . . 

f.  Evapotranspiration.  The  combined  process  of  evaporation  from  the  surface  and 
transpiration  of  water  through  plant  leaves.  The  term  is  also  used  to  denote  the 
amount  of  water  loss  from  the  soil  to  the  air  by  the  evapotranspi  rat  i  on  process. 
E vapot ransp i ra t i on  depends  on  the  evaporative  demand  of  the  air,  moisture 
availability,  soil  type,  and  vegetation. 

g.  Potential  Evapotranspiration.  The  amount  ot  water  lost  from  the  soil  by 
eva po t ra nsp i rat i on  from  the  surface  completely  covered  by  vegetation  with  an  abundant 
water  supply  at  all  times.  Tnornthwaite  (1944)  developed  the  concept  of  potential 
evapotranspiration  to  provide  a  measure  of  the  water  need  of  a  region.  besides 
complete  crop  cover  and  soil  mosture  at  field  capacity,  "the  size  of  the  area  under 
the  high  moisture  conditions  has  to  be  large  enough,  so  that  evapotranspiration  from 
the  area  is  not  affected  by  external  factors  such  as  the  advection  of  moist  or  dry 
air  masses  and  their  modification  by  local  conditions"  ( Thornthwa i te  and  Mather, 
1 9  b  b  ) .  Potential  evapotranspi ration  is  the  maximum  rate  of  evapotranspiration  given 
the  meteorological  conditions  and  is  a  function  of  meteorological  elements  only. 

h.  Moisture  Deficit.  The  difference  between  actual  evapotranspiration  and 
potential  evapotranspiration.  A  positive  moisture  deficit  indicates  the  amount  by 
which  soil  moisture  is  unable  to  support  evapot ra nsp i rat i on  at  the  potential  rate. 


LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


AF  GW  C 

AWN 

E  TP 

NM 

PET 

GPF 

R/E 

USmI l TAG 


Air  Force  Global  Weather  Central 
Automated  Weather  Network 
Evapotransp i rat i on  Potential 
Naut i ca 1  Mile 

Potential  Evapotranspi rat  ion 
Ouantitative  Precipitation  Forecast 
Spread  Ratio 

United  States  Air  force  Environmental  Technical  Applications  Center 
i-Di mens  1 ona I  Nephanalysis  (AFGWC  Automated  Cloud  Analysis  Model) 


30NEPH 


